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Abstract 



We studied the temperature and metal abundance distributions of the intra-cluster medium (ICM) in 
k«J ' a group of galaxies NGC 1550 observed with Suzaku. The NGC 1550 is classified as a fossil group, which 

have few bright member galaxies except for the central galaxy. Thus, such a type of galaxy is important 
to investigate how the metals are enriched to the ICM. With the Suzaku XIS instruments, we directly 
measured not only Si, S, and Fe lines but also O and Mg lines and obtained those abundances to an outer 
region of ~ 0.5 riso for the first time, and confirmed that the metals in the ICM of such a fossil group are 
indeed extending to a large radius. We found steeper gradients for Mg, Si, S, and Fe abundances, while O 
showed almost flat abundance distribution. Abundance ratios of a-elements to Fe were similar to those of 
the other groups and poor clusters. We calculated the number ratio of type II to type la supernovae for 
the ICM enrichment to be 2.9 ±0.5 within 0.1 r 180 , and the value was consistent with those for the other 
groups and poor clusters observed with Suzaku. We also calculated metal mass-to-light ratios (MLRs) for 
Fe, O and Mg with B-band and K-band luminosities of the member galaxies of NGC 1550. The derived 
MLRs were comparable to those of NGC 5044 group in the r < 0.1 rigo region, while those of NGC 1550 
£ — | are slightly higher than those of NGC 5044 in the outer region. 

Key words: galaxies: groups: individual (NGC 1550), galaxies: intergalactic medium, galaxies: 
abundances 



1. Introduction Rcnzini (1997) and Makishima et al. (2001) summa- 
rized iron-mass-to-light ratios (IMLR) with B-band lu- 
Groups and clusters of galaxies play a key role for in- minosity for various objects with ASCA, as a function 
vestigating the formation of the universe and they act of their plasma temperature serving as a measure of the 
as a building blocks in the framework of a hierarchical system richness, and IMLRs in groups were found to be 
formation of structures. The metal abundances of Intra- smaller than those in clusters. They also showed that 
cluster medium (ICM) in groups and clusters carry a lot the early-type galaxies released a large amount of met- 
of information in understanding the chemical history and als which were probably formed through past supernovae 
evolution of groups and clusters. Recent X-ray observa- explosions as shown earlier by Arnaud et al. (1992). In 
tions allow us to measure temperature and metal abun- order to obtain a correct modeling of ICM, we need to 
dance distributions in the ICM based on the spatially re- know the correct temperature and metal abundance pro- 
solved spectra. A large amount of metals of the ICM are files without biases (e.g., Buote 2000; Sanders & Fabian 
mainly produced by supernovae (SNe) in early- type galax- 2002). Especially for the ICM of cooler systems, such as 
ies (Arnaud et al. 1992: Renzini et al. 1993), which are elliptical galaxies and groups of galaxies, careful analysis 
classified roughly as type la (SNe la) and type II (SNe II). are required as mentioned in Arimoto et al. (1997) and 
Because Si and Fe are both synthesized in SNe la and II, Matsushita et al. (2000). 

we need to know O and Mg abundances, which are syn- The spatial distribution and elemental abundance pat- 

thesized predominantly in SNe II, in resolving the past tern of the ICM metals were determined with the large 

metal enrichment process in ICM by supernovae. In order effective area of XMM-Newton (Matsushita et al. 2007b; 

to know how the ICM has been enriched, we need to mea- Tamura et al. 2004; Bohringer et al. 2005; O'Sullivan et al. 

sure the amount and distribution of all the metals from O 2005; Sanders & Fabian 2006; de Plaa et al. 2006; de Plaa 

to Fe in the ICM. et al. 2007; Werner et al. 2006; Simionescu et al. 2008). 
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Table 1. Suzaku Observation logs for NGC 1550. 



Region 


Sequence No. 


Observation date 


(RA, Dec)* 


Exp. 


After screening 








J2000.0 


ksec 


(BI/FI) ksec 


center 


803017010 


2008-08-16T04:27:05 


(04 h 19 m 47?7, +02°24'38") 


83.3 


82.0/83.3 


offset 


803018010 


2008-08-15T02:44:04 


(04 h 2Q m 59?5, +02°24'31") 


41.1 


40.5/40.7 



* Average pointing direction of the XIS, written in the RA_NOM and DECJMOM keywords of the event FITS files. 



On the other hand, the abundance measurements of O and 
Mg with XMM-Newton were possible only for the central 
regions of brightest cooling core clusters due to the rela- 
tively high intrinsic background. Rasmussen & Ponman 
(2007); Rasmussen & Ponman (2009) report the Si and Fe 
profiles of 15 groups of galaxies observed with Chandra. 
They suggest that the Si to Fe ratios in groups tend to 
increase with radius, and the IMLRs within rsoo show a 
positive correlation with total group mass (temperature). 
Suzaku XIS can measure all the main elements from O to 
Fe, because it realizes lower background level and higher 
spectral sensitivity, especially below 1 keV (Koyama et 
al. 2007). Suzaku observations have shown the abun- 
dance profiles of O, Mg, Si, S, and Fe to the outer re- 
gions with good precision for several clusters (Matsushita 
et al. 2007a; Sato et al. 2007a; Sato et al. 2008; Sato et 
al. 2009a; Sato et al. 2009b; Tokoi et al. 2008; Komiyama 
et al. 2009). Combining the Suzaku results with super- 
nova nucleosynthesis models, Sato et al. (2007b) showed 
the number ratios of SNe II to la to be 3.5. 

NGC 1550 is a SO galaxy and one of the nearest (z = 
0.0124) X-ray bright galaxies. The NGC 1550 is also 
classified as a fossil group (Jones et al. 2003), and an X- 
ray extended object RX J0419+0225 was first discovered 
by the ROSAT ALL SKY SURVEY from a position cen- 
tered on the NGC 1550 galaxy. From ASCA observation 
(Kawaharada et al. 2003; Fukazawaet al. 2004), the MLRs 
with B-band is comparable to those of clusters. Sun et al. 
(2003) reports the temperature drop at the central region 
and also declines beyond 0.1 times of the virial radius, 
riso, with Chandra observation. Kawaharada et al. (2009) 
shows the gas mass and metal mass from the temperature 
and metal abundances observed with XMM-Newton. In 
addition, they derived the mass-to-light ratios of O, Si, 
and Fe with near infrared (K-band) luminosity. The re- 
sultant IMLR within ~ 200 I172 kpc exhibits about 2 orders 
of magnitude decrease toward the center. NGC 1550, such 
a fossil group, is a important object to investigate how the 
metals have been enriched to the ICM, because of little 
metal supply from the present-day member galaxies. 

This paper reports on results from Suzaku observa- 
tions of NGC 1550 out to 30' ~ 457 /i^ 1 kpc, correspond- 
ing to ~ 0.47 r 180 . We use H = 70 km s" 1 Mpc" 1 , 
51a = 1 — — 0.73 in this paper. At a redshift of 
z = 0.0124, 1' corresponds to 15.2 kpc, and the virial 
radius, r 180 = 1.95 h^ v /fc(T)/10 keV Mpc (Markevitch 
et al. 1998), is 0.97 Mpc (63f5) for an average temper- 
ature of k(T) = 1.2 keV. Throughout this paper we 
adopt the Galactic hydrogen column density of Ah = 
1.15 x 10 21 cm" 2 (Dickey & Lockman 1990) in the di- 
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Fig. 1. Combined XIS image in the 0.5-2.0 keV energy 
range. The observed XIS0, 1, and 3 images were added on the 
sky coordinate after removing each calibration source region, 
and smoothed with a = 16 pixel ~ 17" Gaussian. Estimated 
components of extragalactic X-ray background (CXB) and in- 
strumental background (NXB) were subtracted, and the ex- 
posure was corrected, though vignetting was not corrected. 
The white circles show the extracted regions in the spectral 
fits. 

rection of NGC 1550. Unless noted otherwise, the solar 
abundance table is given by Anders & Grevesse (1989), 
and the errors are in the 90% confidence region for a sin- 
gle interesting parameter. 

2. Observations and Data Reduction 

2.1. Observations 

Suzaku observed the central and offset regions of 
NGC 1550 in August 2008 (PI: K. Sato). The observation 
logs are given in table 1, and the XIS image in 0.5-2 keV is 
shown in figure 1. We analyzed only the XIS data in this 
paper, because the temperature of the HXD PIN is slightly 
higher during the observations. The XIS instrument con- 
sists of three sets of X-ray CCDs (XIS 0, 1, and 3). XIS 1 
is a back-illuminated (BI) sensor, while XIS and 3 are 
front-illuminated (FI). The instrument was operated in 
the normal clocking mode (8 s exposure per frame), with 
the standard 5 x 5 or 3 x 3 editing mode. During these ob- 
servations, the significant effect of the Solar Wind Charge 
eXchange (SWCX) was not confirmed in ACE data 1 , al- 
though it was known that the SWCX affected the Suzaku 
spectra in the lower energy range as reported in Fujimoto 
et al. (2007) and Yoshino et al. (2009). 

1 http://www.srl.caltech.edu/ACE/ASC 
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Fig. 2. The panels show the observed spectra for the annular regions of NGC 1550 
which are denoted in the panels, and the data are plotted by red and black crosses 
for BI and FI, respectively. The estimated NXB components are subtracted, and the 
black and red lines show the best-fit model for the BI and FI spectra, respectively. The 
BI spectra of the ICM component are shown in orange and magenta lines. The CXB 
components arc shown in a yellow line, and the Galactic emissions are indicated by 
blue and cyan lines, respectively. The energy range around the Si K-edge (1.825-1.840 
keV) is ignored in the spectral fits. The lower panels show the fit residuals in units 
of a. 



2.2. Data Reduction 

We used version 2.2 processing data, and the analysis 
was performed with HEAsoft version 6.6.3 and XSPEC 
12.5.0ac. Here we give just a brief description of the data 
reduction. The light curve of each sensor in the 0.3- 
10 keV range with a 16 s time bin was also examined 
in order to exclude periods with anomalous event rates 
which were greater or less than ±3<r around the mean to 
remove the charge exchange contamination Fujimoto et al. 
(2007), while Suzaku data was little affected by the soft 
proton flare compared with XMM data. The exposure af- 
ter the screening was essentially the same as that before 
screening in table 1, which indicated that the non X-ray 
background (NXB) was almost stable during the obser- 
vation. Event screening with cut-off rigidity (COR) was 
not performed in our data. In order to subtract the NXB 
and the extra-galactic cosmic X-ray background (CXB), 
we employed the dark Earth database by the "xisnxbgen" 
Ftools task. 

We generated two Ancillary Response Files (ARFs) for 
the spectrum of each annular sky region, A v and A B , 
which respectively assumed uniform sky emission and 
~ 1° x 1° size of the /3-model surface brightness profile, 
[3 = 0.47 and r c = 0.'85, in Fukazawa et al. (2004), by 
the "xissimarfgen" Ftools task (Ishisaki et al. 2007). We 
also included the effect of the contaminations on the op- 
tical blocking filter of the XISs in the ARFs. Since the 
energy resolution also slowly degraded after the launch, 
due to radiation damage, this effect was included in the 
Redistribution Matrix File (RMF) by the "xisrmfgen" 
Ftools task. 



3. Temperature and Abundance Profiles 

3.1. Spectral Fit 

We extracted spectra from five annular regions of 0'-3', 
3'-6', and 6'-12' for the central observation, and whole 
area, which corresponds to 12'-30', for the offset region, 
centered on (RA, Dec) = (4 h 19 m 37?9, +02°24'36"). Each 
annular spectrum is shown in figure 2. The ionized Mg, Si, 
S, Fe lines are clearly seen in each region. The O VII and 
O VIII lines were prominent in the outer rings, however, 
most of the O VII line was considered to come from the 
local Galactic emission, and we dealt with those in the 
same way as mentioned in Sato et al. (2007a); Sato et al. 
(2008); Sato et al. (2009a); Sato et al. (2009b). 

The spectra with BI and FI for all regions were fitted 
simultaneously in the energy range, 0.4-7.1/0.4-5.0 keV 
(BI) and 0.6-7.1/0.6-5.0 (FI) for the central/offset obser- 
vations, respectively. In the simultaneous fit, the common 
Galactic emission and CXB components were included for 
all regions. We excluded the narrow energy band around 
the Si K-edge (1.825-1.840 keV) because its response was 
not modeled correctly The energy range below 0.4 keV 
was also excluded because the C edge (0.284 keV) seen in 
the BI spectra could not be reproduced well in our data. 
The range above 7.1 keV was also ignored because Ni line 
(~ 7.5 keV) in the background left a spurious feature after 
the NXB subtraction at large radii. In the simultaneous 
fits of BI and FI data, only the normalization parameter 
was allowed to take different values between them. 

It is important to estimate the Galactic component pre- 
cisely, because the Galactic component gives significant 
contribution especially in the outer regions, as shown in 
figure 2. However, the ICM component is still dominant 
in almost all the energy range except for the OVII line. 
We assumed two temperature apec model (assuming 1 so- 
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3. (a): Radial temperature profiles derived from the spectral fits for each annulus against the projected radius. Black dashed 
show systematic change of the best-fit values by varying the thickness of the OBF contaminant by ±10%. Light-gray dashed 
denote those when the NXB levels are varied by ±10%. (b)-(f): Radial abundance profiles are plotted in the same way as in 



lar abundance and zero redshift) for the Galactic compo- 
nent, and fitted the data with the following model for- 
mula: constant x (apec cool + phabs x (apec hot +ICM)) as 
shown in Yoshino ct al. (2009). As a result, the best- 
fit temperatures of 0.098 and 0.219 keV for the Galactic 
models as shown in table 2 are consistent with the values 
as shown in Yoshino et al. (2009) and Lumb et al. (2002). 
Thus, we concluded that the two temperature model of 
the Galactic emission was enough to represent NGC 1550 
data. The resultant normalizations of the apec models in 
table 2 are scaled so that they give the surface brightness 
in unit solid angle of arcmin 2 , and are constrained to give 
the same surface brightness and the same temperature for 
the simultaneous fits of all annuli 2 . 

The ICM spectra for the central observation, r < 12', 
were clearly better represented by two vapec models than 
one vapec model in the \ 2 test. On the other hand, the 
ICM spectra for the offset observation, r = 12-30', were 
well-presented by a single temperature model. Thus, we 
carried out the simultaneous fit with the following for- 
mula of the Galactic and ICM components: constant x 
(apec cool + phabs x (apec hot + vapec 0<r<30 , + vapec 0<r<12 ,). 
The fit results are shown in table 2. The abundances were 
linked in the following way, Mg=Al, S=Ar=Ca, Fe=Ni, 
which gave good constraint especially for the offset re- 
gions. The abundances were also linked between the two 
vapec components for r < 12' region. Results of the spec- 
tral fit for individual annuli are summarized in table 2 

2 The Normalizations of "normi" in table 5 in Sato et al. (2008) 
are an order of magnitude larger than the actual values. 



and figure 3, in which systematic error due to the OBF 
contamination and NXB estimation are shown. 

3.2. Temperature Profile 

Radial temperature profile and the ratio of the vapec 
normalizations between the hot and cool ICM components 
are shown in figure 3(a) and table 2. The ICM tempera- 
ture of hot and cool components at the central region was 
<~ 1.5 and ~ 0.9 keV, respectively, and the temperature de- 
creased mildly to ~ 1.0 keV in the outermost region, while 
the cool components were almost constant at ~ 0.9 keV. 
Our results for the two temperature ICM model are con- 
sistent with the XMM result (Kawaharada et al. 2009). 
For the hot component, our results are also consistent 
with the previous Chandra result (Sun et al. 2003). The 
radius of 30' ~ 457 kpc corresponds to ~ 0.47 riso, and the 
temperature decline is clearly recognized in such a small 
system out to this radius. 

3.3. Abundance Profiles 

Metal abundances are determined for the six element 
groups individually as shown in figures 3(b)~(f). The four 
abundance values for Mg, Si, S, and Fe and their radial 
variation look similar to each other. The central abun- 
dances lie around ~ 0.5-0.6 solar, and they commonly 
decline to about 1 /4 of the central value in the outermost 
annulus. On the other hand, the O profile looks flatter 
compared with the other elements. Because the results 
for the offset regions had large errors, we examined the 
summed spectra for these regions. We noted that, when 
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Table 2. Summary of the parameters of the fits to each annular spectrum of NGC 1550. All annuli were simultaneously fitted. 
Errors are 90% confidence range of statistical errors, and do not include systematic errors. The solar abundance ratio of angr was 
assumed. These results are plotted in figure 3. 



Galactic 


Normi * 


(keV) 


Norm 2 * 


kT 2 
(keV) 








3.37™ 


0.098iffi 


0-92^;^ 


u.ziy_ 046 






ICM 


l\fn rm i T 


kTi 
(keV) 


1VL/I 1 1 Lj, 


iv-*- 2 

(keV) 


]\fnrm-i I Afnrm^, 

1VL/I 1 1 L\ j 1VL/I 111,2 




0'-3' 

3'-6' 

&-12' 

12'-30' 

total 


358.5t& u 
88.2±f;| 

— 9 
r Q + l.i 


i ' OU -0.02 

1 43+0.03 

i 07+O.O4 
1-0 ' —0.03 

07+0.03 

u ' y ' -0.05 


58.8+t? 

c 4+1.3 
°- 4 -0.9 
r 9 +0.5 

"•^ — 2.8 


0.91^ 

o.86l8:8 

81+°-^ 

u u -0.05 


6.iol8:a 

16.48l»£ 
5 42+2-g 

^ — 2.96 


613/463 
531/469 
520/469 

73/91 
1737/1492 


ICM 




(solar) 


Nc 
(solar) 


Mg,Al 
(solar) 


Si 
(solar) 


S,Ar,Ca 
(solar) 


Fc,Ni 
(solar) 


0'-3' 
3'-6' 
6'-9' 
9'-17' 


2Q+ ' 07 

u - zy -0.07 

n or+0.09 

U ' OO -0.07 

25+ 012 

u - zo -0.09 

00+ ' 47 
u - uu -o.oo 


o.sol^ 

0.26i^ 
0.39l^ 4 


S3+ (U)6 

u ' ol -0.05 
QQ+0-08 

10+ 023 
u - lu -o.io 


S5+ (U)4 

u - oo -0.04 

-5fi+ 04 

U - JO -0.03 
n 9 fi +0.05 
U - zS -0.04 

n+ - 20 

U - 1J -0.13 


^7+ () - ,)4 

u -°'-0.04 

34+ 04 

u -°^-0.04 

31+ 06 

u - ol -0.06 

00+ 25 
u - uu -o.oo 


n 47+0-02^ 
u - 4 ' -0.02 
n 09+0.02 

u - oz -0.02 

29+ 02 

u - zy -0.02 

20+ ' 04 

U ' ZU -0.04 



* Normalization of the apec component divided by the solid angle, Q u , assumed in the uniform-sky ARF 
calculation (20' radius), Norm = J n c n H dV / (4?r (1 + 2) 2 D|) / f2» xlO" 20 cm 5 arcmin 2 , where Da is the 
angular distance to the source. 

t Normalization of the vapec component scaled with a factor of the selected region comparing to the assumed 
image in "xissimarfgen" , Norm = factor J n e n H dV /[47r(l + z) 2 D 2 ] xl0~ 20 cm 5 arcmin 2 , where Da is the 
angular distance to the source, 
t All regions were fitted simultaneously. 



Table 3. List of x 2 /dof for the fits of the nominal and considering the systematic errors such as contaminant of OBF and background 
level. For details, see text. 



Region 


nominal 


contaminant 


background 






+ 10% -10% 


+ 10% -10% 


All 


.. 1737/1492 


1741/1492 1803/1492 


1745/1492 1747/1492 



4. Discussion 

4-1. Metallicity Distribution in the ICM 

Suzaku observation of NGC 1550 confirmed that the 
metals in the ICM of this fossil group are indeed extending 
to a large radius. The measured elements are O, Mg, Si, S, 
and Fe out to a radius of 30' ~ 457 kpc, which corresponds 
to ~ 0.47 riso, as shown in figure 3. The Ne abundance 
had a large uncertainty due to the strong coupling with 
Fe-L lines. Distributions of Mg, Si, S, and Fe are similar 
to each other, while O profile shows no central peak and 
large error in the outer region at r > 12'. We plotted 
abundance ratios of O, Mg, Si, and S to Fe as a function of 
the projected radius in figure 4. The ratios Mg/Fe, Si/Fe 
and S/Fe are consistent to be a constant value around 
1.5-2, while O/Fe ratio in the innermost region (r < 3') 
is significantly lower around 0.5. In addition, the O/Fe 
ratio suggests some increase with radius. Note that these 
abundance profiles are not deconvolved and are averaged 
over the line of sight. 

Recent Suzaku observations have presented abundance 



we examined all abundances to be free in the fits, the re- 
sultant parameters did not change within the statistical 
errors. In addition, even if all regions were fitted by a two 
temperature model, the resultant abundance profiles did 
not change within the statistical errors. 

We also examined the systematic error of our results 
by changing the background normalization by ±10%, and 
the error range is plotted with light-gray dashed lines in 
figure 3. The systematic error due to the background 
estimation is almost negligible. The other systematic error 
concerning the uncertainty in the OBF contaminant is 
shown by black dashed lines as shown in figure 3. A list 
of % 2 /dof by changing the systematic errors is presented 
in table 3. We note that Ne abundance is not reliably 
determined due to an overlap with the strong and complex 
Fc-L line emissions, however we left these abundance to 
vary freely during the spectral fit. 
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Fig. 4. Comparison of the Fe abundance and the other met- 
als to the Fe abundance ratios for NGC 1550 (red crosses) with 
those for the NGC 5044 group (black crosses; Komiyama et 
al. 2009), A 262 cluster (light gray crosses; Sato et al. 2009b) 
and AWM 7 cluster (light gray diamonds; Sato et al. 2008). 
The Fe abundance and the O/Fe, Mg/Fe, Si/Fe, and S/Fc 
abundance ratios in solar units (Anders & Grevesse 1989) are 
plotted against the projected radius scaled by the virial ra- 
dius, ri8o, in all panels. 

profiles in several relaxed cooling core groups and clusters: 
a group of galaxies NGC 5044 (Komiyama et al. 2009), a 
poor cluster of galaxies Abell 262 (Sato et al. 2009b), and 
AWM 7 (Sato et al. 2008). We compare Fe abundance and 
O, Mg, Si, S to Fe abundance ratios of NGC 1550 with 
those of NGC 5044, Abell 262 and AWM 7 as shown in 
figure 4. While NGC 1550 shows slightly lower Fe abun- 
dance than those of NGC 5044, Abell 262, and AWM 7 
in the central region (r ^ 0.1 ?T8o), the abundance in the 
outer region (r > 0.3 riso) are quite similar, showing a de- 
crease to ~ 0.2 solar. On the other hand, the abundance 
ratios of O/Fe, Mg/Fe, Si/Fe, and S/Fe are quite similar 
between the four systems. 

Matsushita (2010) reports the Fe radial profiles of 28 
clusters of galaxies observed with XMM. The Fe abun- 
dance profile of NGC 1550 also has similar feature to 
those of clusters up to ~ 0.5 rigo- Although Rasmussen 
& Ponman (2009) suggests that the radial Si profiles in 
group have softer than the Fe profiles, our results show 
the constant values of the Si to Fe ratios up to ~ 0.5 rigo 
as shown in figure 4. 

4.2. Number Ratio of SNe II to SNe la 

In order to examine relative contributions from SNe la 
and SNe II to the ICM metals, the elemental abundance 
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Fig. 5. Fit results of number ratios of elements to Fe of 
NGC 1550. Top and second panel show the abundance num- 
ber ratios in solar unit in 0.1 < r < 0.2 rigo (black) and0.1ri8o 
(red) regions, respectively. Blue dashed and solid lines in 
the second panel correspond to the contributions of SNe la 
(W7) and SNe II within 0.1 riso, respectively. Ne (atomic 
number = 10) is excluded in the fit. Third panel shows the 
comparison with solar abundance table of Anders & Grevesse 
(1989) (light-gray line) and Lodders (2003) (dashed light-gray 
line). Forth panel indicates ratios of data points to the best-fit 
model. Bottom panel indicates fractions of the SNe la contri- 
bution to total metals in the best- fit model for each element, 
respectively. 

pattern of O, Mg, Si, S and Fe was examined for the in- 
ner (r < 0.1 riso) and the immediate outer (0.1-0.2 r^o) 
regions. The abundance ratios to Fe were fitted by a com- 
bination of average SNe la and SNe II yields per super- 
nova, as shown in figure 5. The fit parameters were the 
integrated number of SNe la (iVr a ) and the number ratio 
of SNe II to SNe la (N u /Ni a ), because iV Ia could be well 
constrained by the relatively small errors in the Fe abun- 
dance. The SNe la and II yields were taken from Iwamoto 
et al. (1999) and Nomoto et al. (2006), respectively. We 
assumed the Salpeter IMF for stellar masses from 10 to 50 
Mq with the progenitor metallicity of Z = 0.02 for SNe II, 
and W7, WDD1 or WDD2 models for SNe la. Table 4 and 
figure 5 summarize the fit results. The number ratios were 
better represented by the W7 SNe la yield model than by 
WDD1. The number ratio of SNe II to SNe la with W7 is 
^2.9 within 0.1 riso, while the ratio assuming WDD1 is 
~ 2.1. The WDD2 model gave the result very similar to 
the W7 value. The resultant number ratios are consistent 
with the previous result by Sato et al. (2007b). We also 
compared the abundance pattern of NGC 1550 with the 
solar abundance. The third panel in figure 5 shows this 
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Table 4. Integrated number of SNc I (AT Ia ) and number ratio of SNc II to SNe la (Ah/Ma)- 
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Region SNc la Model N u N u /Ni a g/dof 



<0.1ri 80 


W7 


1.3 ±0.1 x 10 8 


2.9±0.5 


3.7/3 


0.1 - 0.2 r 180 


W7 


2.6 ± 0.4 x 10 8 


2.8±1.0 


5.7/3 


<0.1r 180 


WDD1 


1.3 ±0.1 x 10 8 


2.1±0.5 


44.2/3 


0.1 - 0.2 r 180 


WDD1 


2.8 ±0.5 x 10 8 


2.0±1.0 


21.0/3 


<0.1ri 80 


WDD2 


1.2 ±0.1 x 10 8 


2.8 ±0.6 


10.4/3 


0.1 - 0.2 r 180 


WDD2 


2.4 ± 0.4 x 10 8 


2.8 ±1.1 


8.9/3 



comparison for r < 0.1 rigo of NGC 1550 with two differ- 
ent solar abundance patterns given by Anders & Grevesse 
(1989) and Lodders (2003). Abundances of Mg, Si, and S 
fall between the two solar abundance patterns. 

Almost ~ 80% of Fe and ~ 40% of Si and S were syn- 
thesized by SNe la in the W7 model, as demonstrated 
in the bottom panel of figure 5. These observed features 
of the fossil group are similar to those for clusters with 
kT = 2 - 4 keV clusters studied by Sato et al. (2007b) and 
Sato ct al. (2009b). The values in table 4 imply that the 
A^ii/iVi a ratio for the inner and outer regions behave in the 
similar manner for different supernova models. We note 
that the fit was not formally acceptable based on the \ 2 
value in table 4. As described in Sato et al. (2007b), the 
models adapted here (SNe yield, Salpeter IMF, etc.) are 
probably too simplified. 

4-3. Metal Mass-to-Light Ratio 

We derived 3-dimensional gas mass profile by extend- 
ing the previous XMM-Ncwton result (Kawaharada et 
al. 2009) for the region within 14' arcmin (<~ 0.2 ri 80 ). 
Combining it with the abundance profiles obtained with 
Suzaku, we calculated cumulative metal mass. The de- 
rived masses of Fe and Mg within the 3-dimcnsional radius 
of r < 457 kpc (r - 0.5 n 80 ) are 1.3 x 10 9 , 3.4 x 10 s M Q , re- 
spectively, and the O mass within r < 183 kpc is 1.7 x 10 9 
Mq. Errors of the metal mass, which were used to calcu- 
late mass-to-light rations, were taken from the statistical 
errors of each elemental abundance in the spectral fits, 
because these are much larger than the error of gas mass 
by Kawaharada et al. (2009). 

We examined mass-to- light ratios for O, Fe, and Mg 
(OMLR, IMLR, and MMLR, respectively) which enabled 
us to compare the ICM metal distribution with the stel- 
lar mass profile. Historically, B-band luminosity has been 
used for the estimation of the stellar mass (Makishima 
et al. 2001), however we calculated it using the K-band 
luminosity in NGC 1550 based on the Two Micron All 
Sky Survey (2MASS) catalogue 3 . This method is useful 
in performing a uniform comparison with the properties 
in other groups and clusters based on the same K-band 
galaxy catalogue to trace the distribution of member el- 
liptical galaxies. 

In the 2MASS catalog, we used all the data in a 2° x 2° 
region around NGC 1550 without the selection of galaxy 
morphology, and subtracted the luminosity in a r > 1° re- 

3 The database address: http://www.ipac.caltech.edu/2mass/ 



gion, which corresponds to about ?T 8 o, as the background. 
We then deprojected the luminosity profile as a function 
of radius assuming a spherical symmetry. In order to con- 
vert the K-band magnitude of each galaxy to the B-band 
value, we assumed the luminosity distance Z?l = 53.6 Mpc, 
and an appropriate color B — K = 4.2 for early- type galax- 
ies given by Lin & Mohr (2004), along with the Galactic 
extinction A B = 0.583 from NASA/IPAC Extragalactic 
Database (NED) in the direction of NGC 1550. 

The integrated values of OMLR, IMLR, and MMLR us- 
ing the estimated B-band luminosity within r < 183 kpc 
(r < 0.2 nso) turned out to be - 3.3 x 10~ 2 , - 7.4 x 10~ 3 , 
and ~ 3.3 x 10~ 3 M q /Lq q , respectively, as shown in ta- 
ble 5. The errors are based only on the statistical errors of 
metal abundance in the spectral fit, and the uncertainties 
in the gas mass profile and the luminosities of member 
galaxies are not included. Note that we did not adjust 
metal-mass and K-band profiles by considering the Suzaku 
PSF effect, because uncertainties in the metal mass had 
the dominant effect in our MLR estimation. 

We compared these B-band MLRs for NGC 1550 with 
those of other groups and clusters. The MLRs are all mea- 
sured within inner (~ 0.1 ri 8 o) and outer (~ 0.25 ?Tso) re- 
gions as shown in table 5 and figure 6. As for Fe (IMLR), 
NGC 1550 shows a similar value with the other groups 
and poor clusters in the inner region. As mentioned in 
subsection 4.1, the Fe abundance itself of NGC 1550 in 
the inner region is slightly smaller than those in other 
groups and clusters. However, the IMLR in this region is 
almost the same as others, due to somewhat low value of 
stellar mass in NGC 1550. Looking at the outer region, 
NGC 1550 shows slightly higher IMLR than NGC 5044 
and comparable to those in other poor clusters. On the 
other hand for the OMLR in the inner region, the poor 
systems show much lower values than the larger high- 
temperature systems. Interestingly, the poor systems (in- 
cluding NGC 1550) show higher OMLR in the outer re- 
gion, comparable to those in the larger systems. The spa- 
tial extent of O looks to be relatively large in these very 
small systems. 

Rasmusscn & Ponman (2009) suggests that the IMLRs 
within rsoo of 15 groups of galaxies observed with Chandra 
have a positive correlation with the groups mass (temper- 
ature). Our results and Chandra's results are also almost 
consistent, and the IMLRs of groups are slightly lower 
than those of clusters. 

We also calculated the MLRs by directly using the K- 
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Table 5. Comparison of IMLR, OMLR and MMLR with B-band luminosity for all systems. 





IMLR 


OMLR 


MMLR 




r 




k(T) 


Reference 




(Mq/L*) 

: ^ll ijj. 


(Mq/L*) 






(kpc) 


(»*i8o) 


(kcV) 




Suzaku 


















NGC 5044 . . . 


2. 6+02 x 10~ 3 


6.6tH x 10~ 3 


1.6i°; 2 x 10- 


-3 


88 


0.10 


- 1.0 


Komiyama et al. (2009) 




3.6±8:| x io- 3 


9A±%i x IO" 3 


2.6lo1 x 10" 


-3 


260 


0.30 






NGC 1550 . . . 


2.6 + ^i x 10~ 3 


1.2l° :3 x IO" 3 


9At\i x 10" 


-4 


46 


0.09 


- 1.2 


This work 




7A+_H x IO" 3 


3.3±o ; g x 10~ 2 


3.3^ x 10" 


-3 


183 


0.19 






Fornax 


4 x 10~ 4 


2 x IO" 3 


- 




130 


0.13 


- 1.3 


Matsushita et al. (2007b) 


NGC 507 ... . 


6.0+H x 10- 4 


2 - 6 -8:l x icr 3 


3.7t° i x 10- 


-4 


120 


0.11 


- 1.5 


Sato et al. (2009a) 




1.7t°* x 10- 3 


6.6±|;| x 10~ 3 


1-lto 2 x 10- 


-3 


260 


0.24 






HCG 62 


2.0±oi x 10~ 3 


6At° i x 10~ 3 


1.0+ ,' 2 x 10- 


-3 


120 


0.11 


- 1.5 


Tokoi et al. (2008) 




4.6+^ x 10~ 3 


3.8±§;I x 10~ 2 


1.5t£$ x 10- 


-3 


230 


0.21 






A 262 


3.6+o ; l x 10- 3 


1.2+° :3 x IO" 2 


1.6 +0 i x 10- 


-3 


130 


0.10 


-2 


Sato et al. (2009b) 




6.7+° 4 x IO" 3 


3.7+\i x 10- 2 


2-7±8:S x io- 


-3 


340 


0.27 






A 1060 


5.7±°r 4 x io- 3 


4.3l°;« x IO" 2 


2.4±°J x 10- 


-3 


180 


0.12 


-3 


Sato et al. (2007a) 




4.0±°,; 4 x IO" 3 


4.3 +2 /° j x IO" 2 


l-6±g;| x 10" 


-3 


380 


0.25 




AWM 7 


4.8+g'l x IO" 3 


2.6 +0 i x IO" 2 


3.4±°j x 10" 


-3 


180 


0.11 


-3.5 


Sato et al. (2008) 




7.6±° A 3 x IO" 3 


x io- 2 


6.7t\i x 10" 


-3 


360 


0.22 




XMM-Newton 


















Centaurus . . . 


4 x 10~ 3 


3 x IO" 2 






190 


0.11 


^4 


Matsushita et al. (2007a) 
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Fig. 6. Comparison of IMLR (a) and OMLR (b) with B-band luminosity to the other clusters and groups within ~ 0.1 riso (black) 
and ~ 0.25 riso (light gray) region. The MLRs of NGC 1550 are shown by the diamonds. 



Table 6. Integrated mass-to-light ratios of O, Mg, and Fc (OMLR, MMLR, IMLR ) with K-band luminosity in units of Mq / Lg . 



Region (kpc/ri 80 ) QMLR MMLR 1 A 1 1 - H 

<45. 7/0.05 9.41 2 - 2 x 10" 4 LlT^ x 10~ 4 2.9j^} x 10" 4 

<91. 3/0.09 3.5±g;£ x IO" 3 2.6±g-J x 10" 4 7.1±g;| x 10" 4 

<182.6/0.19 9.2t° :3 x IO" 3 9.0±^J x 10~ 4 2.0t°, : } x IO" 3 

<456.6/0.47 - 9_j±^5 x 10~ 4 3.5±g;| x 10~ 3 
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Fig. 7. (a) Gas mass, M(< R), and K-band luminosity, Lk within the 3-dimensional radius, R, in NGC 1550, NGC 5044, and 
Abell 262. The gas mass profiles derived from XMM-Newton observations for NGC 1550 (Kawaharada et al. 2009), NGC 5044 
(Komiyama et al. 2009), and Abell 262 (Gastaldcllo ct al. 2007). The red, black, and light-gray lines show each value for NGC 1550, 
NGC 5044, and Abell 262, respectively. The solid and dashed lines corresponds to each gas mass and K-band luminosity, respectively, 
(b) Ratios of the O and Fe mass in units of Mq to the K band luminosity in units of Lq (OMLR and IMLR, respectively) for 
NGC 1550, NGC 5044, and Abell 262, against the 3-dimcnsional radius. The solid and dashed lines corresponds to each OMLR and 
IMLR, respectively. 



band luminosity assuming the Galactic extinction Ak = 
0.050 in the direction of NGC 1550, and the absolute K- 
band solar magnitude of 3.34. The resultant K-band lu- 
minosity within the Suzaku observed region, r < 30', is 
3.8 x 10 11 Lkq, and the radial luminosity profile is also 
plotted in figure 7(a). We calculated the radial profile 
of the OMLR, IMLR, and MMLR values using the K- 
band luminosity out to a radius r ~ 180 kpc (r < 0.2 rigo), 
as shown in figure 7(b) and table 6. The values at the 
outermost radius are ~ 9.2 x 10~ 3 , ~ 2.0 x 10 -3 , and 
~ 9.0 x 10~ 4 Mq/Lkq, respectively. 

In order to investigate the dependence on the system 
size, we compared the MLRs of NGC 1550 (kT w 1.2 keV) 
with those of NGC 5044 (1 keV) and Abell 262 (2 keV), 
which are all poor relaxed cooling-core groups and clus- 
ters. Although they show similar IMLR in the inner re- 
gion, NGC 1550 and Abell 262 show gas mass ratio to the 
K-band luminosity higher than NGC 5044, as shown in 
figure 7(a). As one goes outside (r > 0.1 ri 80 ), NGC 1550 
shows fairly consistent MLRs with Abell 262, larger than 
those in NGC 5044 as shown in figure 6(b). The radial 
profiles of IMLR and OMLR for NGC 1550 look quite 
similar to those for Abell 262 rather than for NGC 5044. 

The radius, r ~ 0.1 rigo, corresponds to the region where 
the gas mass and K-band luminosity seem to overlap as 
shown in figure 7(a). Assuming the stellar mass-to- light 
ratio, Mstar/^K ~ 1 as shown in Arnouts et al. (2007) 
(see also Nagino & Matsushita 2009), this radius indeed 
corresponds to the point where the gas and stellar masses 
are comparable. The similarity of MLRs for the 3 systems 
examined here in the inner region (r < 0.1 rigo) suggests 
that the metal enrichment within the past few Gyr has 
occurred in a similar way. On the other hand, in the 
outer region (r > 0.1 tt 80 ), NGC 5044 shows lower IMLR 
than those in NGC 1550 and Abell 262, even though Fe 
abundances themselves are comparable in these regions. 



This suggests that distribution of stellar mass has some 
different history between these systems. 

We will consider how the observed features of abun- 
dance and MLR profiles can constrain enrichment sce- 
nario. First, metals in the inner region (r < 0.1 rigo) have 
been mostly supplied by the central galaxy, so this re- 
gion should be set aside in the present discussion. In the 
immediate outer region (r > 0.1 rigo), metals from the 
central galaxy could not reach in a few Gyr time scale. 
Also, we may assume that almost the same amount of 
metals per stellar mass was synthesized in all systems be- 
fore the collapse of groups and clusters. In this case, we 
expect very similar MLRs in different systems, contrary 
to the observed feature. As shown in figure 7(b), at least 
NGC 5044 shows lower IMLR and OMLR profiles com- 
pared with those for two other systems. This implies that 
the thermal and/or dynamical evolutions of the gas have 
different history among different systems during or after 
the collapsing period of individual groups and clusters. 

Renzini (2005) showed the expected MLRs as a func- 
tion of the IMF slope. As for OMLR, the expected value 
to be ^0.1 Mq/Lb at a Salpeter IMF is slightly higher 
than our results of ~ 0.03 Mq/Lb within ~ 0.25 rigo as 
shown in figure 6(b). However, as shown in figure 7(b), 
because the OMLRs look increasing toward outer region, 
the OMLRs within the virial radius would be represented 
with the Salpeter IMF. Renzini (2005) also suggested that 
a top heavy IMF would overproduce metals by more than 
a factor of 20, which is much larger than the observed 
values including our results. 

We stress that high-sensitivity abundance observation 
to the outer region of clusters will give important clues 
about their evolution. If O distribution, as well as Fe, 
could be measured to the very outer region (r ~ rigo), 
we may obtain a clear view about when O and Fe were 
supplied to the inter galactic space because most of O 
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should have been synthesized by SNe II and supplied in 
the starburst era. Another possibility is very early metal 
enrichment of O by galaxies or massive Population III 
stars (e.g.Matteucci et al. 2006) before groups and clus- 
ters assemble. In this case, a large part of the intergalac- 
tic space would be enriched quite uniformly with O and 
other elements. Mctallicity information in cluster out- 
skirts would thus give us unique information about the 
enrichment history. For this purpose, instruments with 
much higher energy resolution, such as microcalorimeters, 
and optics with larger effective area will play a key role in 
carrying out these studies. 

5. Summary and conclusion 

Suzaku observation of the fossil group NGC 1550 
showed spatial distributions of temperature and metal 
abundances for O, Mg, Si, S, and Fe up to ~ 0.5 ri 80 
for the first time, and confirmed that the metals in the 
ICM of this fossil group are indeed extending to a large 
radius. The ICM temperature decreases mildly from ~ 1.5 
keV to ~ 1.0 keV in the outer region, similar to the feature 
seen in other clusters. The abundances of Mg, Si, S, and 
Fe drop from subsolar levels at the center to ~ 1/4 solar 
in the outermost region, while the O abundance shows a 
flatter distribution around ~ 0.5 solar without the strong 
concentration in the center. The abundance ratios, O/Fe, 
Mg/Fe, Si/Fe, and S/Fc for NGC 1550 are generally sim- 
ilar to those in groups and poor clusters. The abundance 
pattern from O to Fe enabled us to constrain number ra- 
tio of SNe II to la as 2.9 ±0.5, which is consistent with 
the values obtained for other groups and clusters. The de- 
rived MLRs of NGC 1550 using the B-band and K-band 
luminosities are consistent with those in the NGC 5044 
group for the inner region r ^ 0.1 riso, while NGC 1550 
shows slightly higher MLRs than NGC 5044 in the outer 
region. This suggests that metal enrichment process may 
reflect the size of the system in the sense that larger sys- 
tems contain higher amount of metals for a given stellar 
mass. 

Authors thank the referee for providing valuable com- 
ments. K.S is supported by a JSPS Postdoctral fel- 
lowship for research abroad. Part of this work was 
financially supported by the Ministry of Education, 
Culture, Sports, Science and Technology, Grant-in- Aid for 
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